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This paper is devoted to consideration of the hard exclusive processes e + e~ — > VP, where V = 
p,4>; P — rj, 77'. Experimental measurement of the cross section of the process e + e~ — > (fir] at BaBar 
collaboration at large center mass energy %/s = 10.6 GeV and some low energy experimental data 
\fs ~ 2 — 4 GeV give us the possibility to study the cross section in the broad energy region. 
As the result, we have determined the asymptotic behavior of the cross section of e + e~ — > (fir] in 
the limit s — > 00, which is in agreement with perturbative QCD prediction. Assuming that the 
same asymptotic behavior is valid for the other processes under consideration and using low energy 
experimental data we have predicted the cross sections of these processes at energies y/s — 3.67, 10.6 
GeV. In addition, we have calculated the cross sections of these processes at the same energies within 
perturbative QCD. Our results are in agreement with available experimental data. 

PACS numbers: 12.38.-t, 12.38.Bx, 13. 66. Be, 



Exclusive hadron production in high energy electron-positron annihilation is a very interesting task for theoretical 
and experimental investigations. The presence of high energy scale -y/s that is much greater than typical hadronic 
scale allows one to separate the amplitude of such processes into hard part (creation of quarks at very small distances) 
and soft part (subsequent hadronization of these quarks into experimentally observed mesons at larger distances). 
The first part of the amplitude can be calculated within perturbative QCD. The second part of the amplitude is 
described by distribution amplitudes (DA), which contain nonperturbative properties of final hadrons. 

The description of hard exclusive hadron production within this pattern gives some very interesting predictions of 
the properties of hard exclusive processes 0,0]. One of such prediction is the asymptotic behavior of the amplitudes 
and cross sections of hard exclusive processes in the limit s — > 00 @, 0, [B| • It turns out that this behavior is determined 
by the perturbative part of amplitude and quantum numbers of final hadrons and does not depend on DAs of final 
hadrons. 

To give quantitative prediction for the cross section of hard exclusive process one needs to know DAs. It is interesting 
to note that if DAs are known the theory can equivalently well predict the cross sections for the production of hadrons 
composed of light (u, d, s quarks) or heavy quarks (b, c quarks) . There is well known example of exclusive process 
with heavy quarkonia production e + e~ — > J/^Vc measured at Belle fo] and BaBar [?] collaborations. This process 
was extensively studied in many papers [1, H, Eft El, E2, E, EH, El, Eft E3] within different approaches, what led to 
a better understanding of charmonia properties and production processes. It is important to note that the approach 
to hard exclusive processes described above leads to a reasonable agreement with the experiments. 

In this paper we study the processes e + e~ — * VP, where V = p, <fi\ P = r],r]' ', Experimental measurement of the 
cross section of the process e + e~ — > (fir] at BaBar collaboration at large center mass energy y/s — 10.6 GeV and 
some low energy experimental data y/s ^2 — 4 GeV allow us to study the cross section of this process in the broad 
energy region. Our first purpose is to use these data in order to determine the asymptotic behavior of the cross 
section of the process e+e~ — > (fir]. Assuming that the same asymptotic behavior is valid for the other processes 
under consideration and using low energy experimental data one can predict the cross sections of the e + e~ — > VP at 
energies y/s — 3.67, 10.6 GeV. In addition, we apply perturbative QCD approach to estimate the values of the cross 
sections for the processes under consideration. 

This paper is organized as follows. In the next section we analyze the experimental data for the process e + e~ — > (^77 
and determine the asymptotic behavior of this cross section. Then we apply the result of this study to predict the cross 
sections of the other processes under consideration at the center mass energies y/s = 3.67, 10.6 GeV. In section Hill we 
give theoretical predictions for the cross sections c(e + e~ — > p°rj), <r(e + e~ cr(e + e~ — > (fir)) and a(e + e~ — > (fir)') 
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Figure 1: Typical diagrams for the e + e — > VP process. 



at y/s = 3.67 GeV and 10.6 GeV and compare them with available experimental data. The final section is devoted to 
the discussion of the results of this paper. 



II. THE ASYMPTOTIC BEHAVIOR. 



The amplitude of the process involved can be written in the following form: 

M(e+e-^VP) = 47 ra g(gl)7Mg2) (Vfr, A)P(p,) | JT\ 0) 



•s 



where a is the electromagnetic coupling constant, u(q 2 ) and v(qx) are electron and positron bispinors, s = (qi + q2) 2 
is the invariant mass of e + e~ system squared and J^ m is electromagnetic current. The matrix element (VP \J^ m \ 0) 
can be parameterized by the only formfactor F(s): 

(V( Pl ,X)P(p 2 )\j; m \0) = ie^ a pelp<?plF(s), (1) 

where t\ is the polarization vector of meson V. The cross section of the process under consideration equals 

ira 2 /2|p| \ 3 , , Sl2 



a(e + e- -> VP) = — I J |F( S )| Z • (2) 

In the last formula p is the momentum of the vector meson V in the center mass frame of final mesons. 

In this section we will be interested in the asymptotic behavior of the formfactor F(s) in the high energy region. 
Typical diagrams of the process under consideration are shown in Fig. [TJ The asymptotic behavior of the diagram 
shown in Fig. QJi is a ~ 1/s 2 . At extremely large energies this diagram gives the dominant contribution. However, 
the amplitude of this diagram is suppressed by the smallness of the electromagnetic coupling constant and our study 
shows that in the energy region analyzed in this paper the contribution of this diagram is negligible. Further let us 
consider the diagrams shown in Fig. [T}3,c. According to perturbative QCD |3| amplitudes {1} from such diagrams 
have the following asymptotic behavior 

|Ai + A 2 | + l 

(H 1 { Pl ,X 1 )H 2 ( P2 ,X 2 )\j; m \0) 

where Hi and H 2 are mesons with momenta pi , p 2 and helicities Ai and A2 . For the process under consideration H% 
and H 2 are the vector and pseudoscalar mesons respectively. The helicity of the pseudoscalar meson is, obviously, 
A 2 = 0. Because of antisymmetric tensor in Jl]), longitudinal polarization of the vector meson is forbidden and it is 
transversely polarized (Ai = ±1). So, we have F(s) ~ 1/s 2 , er(e + e~ — > VP) ~ l/.s 4 . On the other hand, in papers 
[Isl. fl9|. it was stated that experimental data can be described only by the dependence a ~ 1/s 3 . To clarify this 
situation we will parameterize the formfactor F(s) by the expression 

F(s) = ^« (3) 
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where a n (s) slightly depends on s due to power and logarithmic corrections to the leading asymptotic behavior. Later 
in this section we will neglect such dependence. This approximation allows us to fix the constants a n from the low 
energy data and predict the cross sections at yfs = 3.67 GeV and 10.6 GeV. We check three different hypothesis 
(n = 3,4 and 5) and compare the results with existing experimental data. 

First we are going to consider the process e + e~ — > cjyq. To fix the constants a n one can use low energy data [20]. 
When the constants a n for different hypothesis are fixed, one can use them to predict the cross section in the high 
energy region and then compare this prediction with available high energy data [2l| measured by BaBar collaboration 

0-BaBar(e + e- -> (f)!]) = 2.9 ± 0.5 fb. 

Our results are shown in Fig. [2]and Table[H From these figure and table it can be clearly seen, that only n — 4 describes 
satisfactory low and high energy experimental BaBar data [13, HH]. This result is in agreement with the predictions 
of perturbative QCD. Form Fig. [2] one can see that low energy CLEO-c data [22] are in disagreement with hypothesis 
n = 3,4. Only hypotheses n = 5 does not contradict to the experimental results. However, if we assume that this 
hypotheses is correct we will faced with dramatic contradiction with the high energy BaBar data (see Fig. [2]and Table 
H]). It should be also noted that in papers [3, Ell it was stated that the energy dependence cr(e + e~ — > c/rq) <~ 1/s 3 
describes experimental data more accurately. We believe that the disagreement of this statement with our conclusion 
arises from the fact that the authors of papers fisi fl9| did not take into account low energy experimental result [20| . 

In view of this the question arises: is it possible to use asymptotic behavior of the cross section ([3]) in the region 
y/s S (2, 3.5) GeV. First, one can estimate the cross uncertainty due to the power corrections as ~ M 2 / s ~. Even, for 
the heaviest meson 4> and the smallest y/s from the reg ion y/s € (2, 3.5) the error is ~ 0.25%. We can also determine 
the size of power corrections from the fitting data [20j by the asymptotic form |(3]) plus next-to-leading-order power 
correction. Our analysis shows that the uncertainty is not greater than 10% ~ 20%. We believe that these arguments 
confirm the applicability of the asymptotic expression for the cross section. 

Now let us consider the processes e + e~ — > pi], e + e~ — ► pry'. Unfortunately at the moment only the low energy 
CLEO-c data [12] are available for these processes. As it is seen from Tab. Uthe experimental error of these data are 
rather large. More precise experimental data can be obtained from the decays J/ip — ► pr], J/ip — > pj]' . Corresponding 
branching fraction are equal to [23] 

Br (e+e" -> pr]) = (1.93 ± 0.23) x 10~ 4 , 
Br (e+er -> pif) = (1.05 ± 0.18) x 10" 4 . 

Generally speaking, these decays can proceed both via strong and electromagnetic interaction (see Fig. [3] for the 
typical diagrams). Because of the isospin violation, however, the gluon induced diagrams are strongly suppressed, and 
purely electromagnetic diagram gives the dominant contribution. The branching fraction of these decays are equal to 

2 

M^Br(J/V; -> e+e"). (4) 

From this relation we can determine the value of the function a n (M 2 ^) and, neglecting energy dependence in a n (s) 

(this assumptions leads to additional errors, which can be estimated as ~ M 2 /s ~ 10%), predict the cross sections 
of the processes e + e~ — > pi] and e + e~ — » ,077' over the large energy region. The results are shown in Figs. IH[5]and 
Table HI From these results it is seen that the low energy data do not allow to understand what hypotheses gives the 
best agreement with the experiments. So, to determine the asymptotic behavior unambiguously the high energy data 
are needed. However, if we assume that the asymptotic behavior for the cross section of the process e + e _ — * (pr) is 
the same as that for the processes e + e~ — * prj, e + e~ — > pr)', one can predict the cross sections of these processes at 
energies yfs = 3.67, 10.6 GeV 

°VJ=3.67 Gev( e+e ~ -> PV) = 8 ± 2 pb, cr^=3.67 Gev( e+e ~ PV ) = 5 ± 3 pb, 

°"v^=io.6 GeV (e + e _ pi]) =2.1±0.5fb, <J^= W . 6 GeV (e + e" -> pr]') = 1.4 ± 0.4 fb. (5) 

The cross sections at the energy y/s = 3.67 GeV are in agreement with the CLEO-c data. 

Let us now proceed with the reaction e + e~ — > <jyq' . In this case experimental information is even more poor. In the 
case one cannot use the branching fraction of the decay Br(J/4 r — > <jyq'), since in this case isospin is conserved and the 
diagram shown in Fig. [3^, gives the main contribution to J/ip — > <\yrj decay. In addition, CLEO-c gives us only upper 
bound on the cross section er(e + e~ — > <f>rf), and no experimental information in high energy region is available. Using 
rj — rf mixing it is possible, however, to estimate the cross section <r(e + e~ — * tfyrf) from the value of the e + e _ — > cfrq 
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Figure 2: Different hypothesis on the energy dependence for the cr(e + e~ — » <frr)). The constants a n are fixed from the experi- 
mental values of this cross section in the low energy region (2 GeV < \fs < 3.5 GeV). 




Figure 3: Typical diagrams for J /if] — > VP decay 
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Figure 4: Different hypothesis on the energy dependence for the a(e + e — > prf). The constants a n are fixed from the J/tp — » pr) 
branching fraction (the rightmost point represents the cross section at y/s = Mj/^ calculated using relations J2}, ((3j) and |[4|). 
) 
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Figure 5: Different hypothesis on the energy dependence for the a(e + e — > prf). The constants On are fixed from the J/ip — > pr/ 
branching fraction (the rightmost point represents the cross section at V* = Afj/,/, calculated using relations J2|), J3} and 
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cross section. The mixing of pseudoscalar mesons can be described by different parameterizations [24]]. In our paper 
we will use the parametrization of rj — rj mixing in quark flavor basis with one mixing angle [H[ : 



i) \ _ / cos <f> — sin <f> \ / t] n 
rj J ~~ I sin $ cos $ J I rj s 



(6) 



where r\ n = (uu + dd)/^/2 and r) s = ss represents the basis of the quark mixing scheme. Using this mixing scheme it 
is easy to obtain the following relation between e + e~ — > (f>rj and e + e~ — > 4>rj' cross sections: 

cr(e + e" — > Sri') 7 
— — — = cot $. 

a(e + e~ — > </>ry) 

There are plenty of theoretical and experimental works dedicated to the determination of the mixing angle. The well 
known estimation, based on Gell-Mann-Okubo mass formulas give the value of $ about 32° for linear mass formula 
and about 45° for a quadratic case. The analysis of the axial anomaly generated decays rj, rf — ► 77 was performed in 
papers [IE H3| and the estimate $ = 30° -j- 35° was obtained. In papers [H, another anomaly based investigation 
of a large set of decay processes was performed and the value <f> = 38° ± 2° was presented. The authors of the recent 
work [30j used dispersive approach to 77 — rj' mixing and obtained the value 

$ = 39.4° ±1°. 



It is interesting to note that this value is close to phenomenological value $ = 39.3° ± 1°, presented in pioneering 
work [25]. We will use this value of the mixing angle in our article. The error of this angle is small in comparison 
with the error in e + e~ — > 4>rj cross section. 

Knowing the value of the cross section a(e + e~ — * cfyq) one can determine the value of the cross section cr(e + e _ — > <j>rf) 
at the energy y/s = 10.6 GeV. 

^=10.6 Gev(e + e- -> 0t/) = 4.2 ± 0.7 fb. 

From this we can calculate the value of the constants a n for different hypothesis and predict the cross section of the 
process e + e~ — > <prf in the low energy region. We show the energy dependence of this cross section in Fig. [6] From 
this figure one sees that the value n = 5 contradicts experimental results at y/s = 3.67 GeV, while n = 3 and n = 4 
do not. From theoretical arguments and presented above analysis of e + e~ — » <prj reaction we think, that the value 
n = 4 is more preferable. For this hypothesis we can predict the value of the cross section at energy yfs = 3.67 GeV 

<V?=3.67 Gev(e + e- -> H) = ^.8 ± 2.1 pb, 

which does not contradict to the CLEO-c data. 

In table]!] we present the results of this section. Second column contains constants 03 4 5 for different final states, 
obtained from low energy fits. In the third and fifth columns experimental results for the cross sections er(e + e~ — * VP) 
at the center mass energies y/s = 3.67 GeV and 10.6 GeV are presented. The results of the calculation are shown in 
the forth and sixth columns. From this table it is clear, that only relation a ~ 1/s 4 agrees with experiment. 



III. CALCULATION OF THE CROSS SECTIONS. 
A. Numerical parameters and distribution amplitudes. 



Now let us try to estimate the cross sections of the processes studied in the last section theoretically. To calculate 
the cross sections of the processes e + e~ — > VP for V — (p, p, P — 77, 77' one needs to know distribution amplitudes 
(DA) of final mesons. For the vector mesons the DAs needed in the calculation can be written as 0, E3] 

(Vm*>V)*>lr*)\% = ^j\xe^*-v{p^V L {x) + (e x -p^j V±(x) + 

l^{a^p»)V T {x) + \{e, vaPl , l5 e\p a zP)V A {x)X , (7) 

where x is the fraction of momentum carried by quark, fv, frit 1 ), Mv are the leptonic, tensor constants and the mass 
of vector meson. In the calculation 4> meson is assumed to be composed of s quarks, so in (7]) q — s. For p° meson 
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ct(a/s = 3.67 GeV) pb 


a(Vs = 10.6 GeV) fb 


VP 


parrams 


exp 


fit results 


exp 


fit results 


m 


a 3 = 1.8±0.2GeV 2 
a 4 = 5.6±0.7GeV 3 
a 5 = 17.3 ± 2.1 GeV 4 


10 ±2.5 


12 ±3 
8± 2 
6± 1 




24 ±6 
2.1 ±0.5 
0.18 ±0.04 


PV 


a 3 = 1.5 ±0.2 GeV 2 
a 4 = 4.7 ± 0.6 GeV 3 
a 5 = 14.5 ± 1.9 GeV 4 


2.1 ± 1.6 


7± 2 
5±3 
3.5 ±0.9 




17 ±4 
1.4 ±0.4 
0.12 ±0.03 




a 3 = 2.7 ± 0.2 GeV 2 
a 4 = 6.4 ±0.4 GeV 3 
a 5 = 14.8 ± 1.1 GeV 4 


2.1 ± 1.9 


23 ±4 
9.8 ± 1.3 
3.8 ±0.5 


2.9 ±0.5 


52 ±8 
2.7 ±0.4 
0.13 ±0.02 


H 


a 3 = 0.76 ± 0.06 GeV 2 
a 4 = 8.1 ±0.7 GeV 3 
a 5 = 85. ± 7. GeV 4 


< 12.6 


1.5 ±0.2 
12.8 ±2.1 
110 ± 18 




4.2 ±0.7 
4.2 ±0.7 
4.2 ±0.7 



Table I: The constants a„ and the cross sections at y/s — 3.67 and 10.6 GeV in comparison with the experimental data. The 
second column contains the constants 03,4,5 for the different final states, obtained from the low energy data. In the third and 
fifth columns the experimental results for the cross sections cr(e + e~ — > VP) at the center mass energy y^s = 3.67 GeV and 
\fs = 10.6 GeV are presented. The results of the calculation are shown in the forth and sixth columns. 

isospin 1 combination qq — (uu — dd)/y2 is assumed. The models for DAs and the decay constants fv, /t(m) f° r 4> 
and p mesons will be taken from paper [31(. 

In the framework of the quark mixing scheme (jGj) the decay constants 

(P(p)\nn 5 n\0) = iftff, 
(P(p)\s^ 75 s\0) = iftpP, 

needed in the calculation can be expressed through the constants 

(Vn(p)\ni^ 5 n\0) = t f n p^ 
( Vs (p)\srj 5S \0) = i] s p^ 

as follows 

% fn \ = ( cos$ ~sin$\ ( P \ 
# ) [ sin$ cos* J [ f J ' 

In turn, the constants /", f s and the mixing angle $ can be determined from experiment [25j 

/" = (1.07 ±0.02)/,, / s = (1.34 ±0.06)/,. 
Within this mixing pattern the DAs needed in the calculation can be written in the following form 0, [l3| 

(P(p)\nf3(z)n a (-z)\0)^ 
(P(p)\sp(z) Sa (-z)\0)^ 

where 



mi) cos $ sin $ , 
to 2 ) cos $ sin <£> . 

The calculation will be done with the masses m s (l GeV) = 150 MeV, m n (lGeV) = (m„(l GeV) + m d (l GeV))/2 = 5 
MeV and it will be assumed that P%{y) = Pa{v) = Pa{u) and P P (y) = Pp{y) = Pp{y)- The models of the leading 



= f d ye^y-v{t^ P2(y) - /» 75 Pp(y)], 

= ^f% iW!SH) {|^(!/)-/;M^p(4 



2m n (fi) 
1 

2m s (//) 



to 2 cos $ + mil sin $ 



to 2 , cos 2 <I> ± ml sin 2 $ 



r 

V2f s 



(to 2 , 
(to 2 , - 
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twist DAs Pa for the r) and rf mesons will be taken from paper [32]. For the function Pp(y) the asymptotic form will 
be used. 

It should be noted that the wave functions and the constants introduced above depend on renormalization scale /x. 
Our calculation shows that the scale dependence of the DAs is not very important and below it will be ignored. At 
the same time the scale dependence of the constants is important and it will be taken into account. The calculation 
will be done at the scale /U = y/s/2. 



B. Numerical results and discussion. 



Having introduced the designations of the DAs one can proceed with the calculation of the cross sections. First we 
are going to consider the diagrams similar to that shown in FigQJi. It is not difficult to calculate the contribution of 
these diagrams to the formfactor F(s) 

\F(s)\ = — + f^PA{y), 
s My V V 2 J Jo 2/12/2 

where e u , e^, e s are the charges of u, d, s quarks correspondingly. For the process e + e~ — > 4>r] at the energy y/s = 10.6 
GeV we have a — 0.01 fb, what is by two orders of magnitude less than the experimental result. Form this it is clear 
that the contribution of the diagrams shown in FigQJi is negligible and below it will be ignored. It should be noted 
that in the limit s — > oo the contribution of FigQJi diagrams to the cross section has the following behavior ~ 1/s 2 . 
So, at extremely large energy this diagrams give the dominant contribution. Our calculation shows that this happens 
at energies y/s > 35 GeV. 

Now we are going to consider the other diagrams shown in Fig. [TJ The leading asymptotic behavior of these 
diagrams is a ~ 1/s 4 . First, it should be noted that the contribution of the diagrams shown in Fig. [TJd is very small 
[l8T | and it will be ignore below. This fact results from rather small admixture of \GG) fock state in the pseudoscalar 
mesons T), rj '. Now we proceed to the calculation of the diagrams shown in Fig. [lfc. The contribution of these diagrams 
to the F(s) for the processes under study can be written as follows [13] 

\F*,(«)\ 
\F W (s)\ 

\F PV (s)\ 



32tt 


f v f*M v sm<f> 


e s Ia, 


9 


s 2 


32tt 


fvf s M v cos$ 


eJo, 


9 


s 2 


32tt 


fvf n M v cos$ 


(e„ - 


9 


2s 2 


32tt 


/y/"AfySin$ 


(e« - 


9 


2s 2 



Fpr,'{s)\ = — r-5 (e«-ed)J . (8) 



In the above expressions 



/*(/*),, ,V T (x)P P (y) 



I Q = / dx dy a s (fi) <——f p (ii) = h (9) 

J a Jo { M v x 2 y 

l V L (x)P A (y) . , V±(x)P A (y) , 1 (l + y)V A (x)P A (y) \ 
2 xy xy 8 x y J 

where ft(fx) — fr(fJ-)/fv, fp{p) — fpin) for the processes with <p meson in the final state and f p {p) = fp{n) for the 
processes with p meson in the final state. 

The following point deserves consideration. The models for the DAs that we use in our calculation are truncated 
series in Gegebauer polynomials. This means that the end point behavior (x — * 0, 1) of these DAs coincides with the 
end point behavior of the corresponding asymptotic DAs. From this it is not difficult to see that the integral Io is 
logarithmically divergent. One way to regularize this divergence is to introduce cut off parameter Xq. In our calculation 
we take the following value of cut off parameter xq = A/y/s. A is of order of typical hadronic scale, which is of order 
of several hundreds MeV. Physical meaning of this cut off can be understood as follows: if x ~ xo quark momentum is 
of order of ~ A and in this region one must take into account transverse motion in hadron what regularizes the whole 
integral Iq- The parameter A can be determined from available experimental results. We determined this parameter 
from the cross section of the process e + e~ — > §r\ at yfs — 10.6 GeV. Thus we get A = 130^ 2 y MeV. The variation of 
parameter A corresponds la deviation from the central value measured at Belle collaboration. 

The results of the calculation are shown in Table HH The second and fifth columns contain the experimental results 
for the cross sections at the energies y/s — 3.67 GeV and y/s = 10.6 GeV correspondingly. The result obtained in 
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a(y/s = 3.67 GeV) pb 


a(Vs = 10.6 GeV) fb 


VP 


exp [22] 


[18] 


this work 


exp [21] 


[18] 


this wok 


m 


10 ±2.5 


8.1-M6.6 


3.7- 


-6.3 




2.4- 


-3.1 


2.4- 


-3.5 


m 


2.1 ± 1.6 


4.3-^8.6 


2.1 - 


-3.6 




1.5- 


- 2.1 


1.6- 


-2.3 




2.1 ± 1.9 


9.6 -f- 19.1 


3.7- 


-6.1 


2.9 ±0.5 


3.3- 


-4.3 


2.4- 


-3.4 


4>v' 


< 12.6 


11.5-^22.6 


4.6- 


-7.6 




4.4- 


-5.8 


3.5- 


-5.0 



Table II: The results of the calculation. The second and fifth columns contain experimental results for the cross sections at 
energies \fs = 3.67 GeV and \fs — 10.6 GeV correspondingly. The result obtained in paper [isl | are shown in the third and 
sixth columns. The results obtained in this paper are shown in the fourth and seventh columns. 

paper [18] are shown in the third and sixth columns. The results obtained in this paper are shown in the fourth and 
seventh columns. The variation of the results are due to the variation in parameter A. It is seen from this table that 
the results of the calculation are in satisfactory agreement with the experiment. 

It should be noted here that in addition to the diagrams shown in Fig. [IJ there is additional contribution to 
the formfactor F(s) which was not considered in this paper. This contribution appears if one takes into account 
higher fock state of the vector and pseudoscalar mesons JogG) and it's asymptotic behavior is also a ~ 1/s 4 . This 
contribution was considered in papers [2. l33j|. In paper [2] the authors asserted that the fock state \qq~G) gives very 
important contribution to the cross section and must be taken into account. Contrary to this conclusion, the author of 
paper [33] asserted that the main contribution arises from the diagrams shown in Fig. Hfc. So, the question about the 
role of higher fock state \qq~G) in the total cross section deserves separate consideration and it will not be considered 
here. 



IV. DISCUSSION AND CONCLUSION. 

In this work we have studied the production of light mesons prj, prj' , (\>r\ and 4>rj' in the high energy electron-positron 
annihilation. 

The first question studied in this paper is the asymptotic behavior of the cross sections of the processes under 
consideration. Perturbative QCD predicts, that the cross section of the reaction e+e~ — * VP has the following 
asymptotic behavior a ~ 1/s 4 in the limit s — > oo. Experimental measurement of the cross section of the process 
e + e~ — > cf>r] at the large center mass energy t/s = 10.6 GeV [2l| and the low energy experimental data y/s ~ 2 — 4 
GeV [20] give us the possibility to study the cross section in the broad energy region. As the result, we have 
determined the asymptotic behavior of the cross section of e + e _ — > (fyq in the limit s — » oo, which is in agreement 
with perturbative QCD prediction. As to the other processes under study, there are no high energy experimental 
data, which allow us to confirm perturbative QCD prediction for these processes. We would like to stress here that 
the high energy experimental data turned out to be crucial in the determination of the asymptotic behaviour of the 
cross sections. Assuming that the asymptotic behavior predicted by perturbative QCD is valid for the other processes 
under consideration, we have calculated the cross sections of the processes e + e~ — > pr),pr)',(f>ri,<fir)' at the energies 
= 3.67, 10.6 GeV. 

In addition, we have applied perturbative QCD approach to calculate the cross sections of the processes under 
study at the energies y/s = 3.67 GeV and y/s — 10.6 GeV. The results of this calculation are in satisfactory agreement 
with available experimental data. 
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